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Abstract 

When interest rate dynamics are described by the Libor Market Model as in Brace, Gatarek & Musiela (1997), 
we show how some essential risk-management results can be obtained from the dual of the calibration program. In 
particular, if the objective is to maximize another swaption's price, we show that the optimal dual variables describe a 
hedging portfolio in the sense of Avellaneda &. Paras (1996). In the general case, the local sensitivity of the covariance 
matrix to all market movement scenarios can be directly computed from the optimal dual solution. We also show how 
semidefinite programming can be used to manage the Gamma exposure of a portfolio. 

Keywords: Libor Market Model, Inverse problems, Semidefinite Programming, Calibration. 

1 Introduction 

A recent stream of works on the Libor Market Model have showed how swap rates can be approximated by a basket 
of lognormal processes under an appropriate choice of forward measure. This, coupled with analytic European basket 
call pricing approximations, allows to cast the problem of calibrating the Libor Market Model to a set of European 
caps and swaptions as a semidefinite program, i.e. a linear program on the cone of positive semidefinite matrices (see 
Nesterov & Nemirovskii (1994) and Boyd & Vandenberghe (2004)). This work exploits the related duality theory to 
provide explicit sensitivity and hedging results based on the optimal solution to the calibration program. 

The lognormal approximation for basket pricing has its origin in electrical engineering as the addition of noise 
in decibels (see for example Schwartz & Yeh (1981)). Its application to basket option pricing dates back to Huynh 
(1994) or Musiela & Rutkowski (1997). Brace, Dun & Barton (1999) tested its empirical validity for swaption pricing 
and Brace & Womersley (2000) used it to study Bermudan swaptions. More recently, d'Aspremont (2003), Kawai 
(2002) and Kurbanmuradov, Sabelfeld & Schoenmakers (2002) obtained additional terms in the expansion and further 
evidence on the approximation performance. On the calibration front, Rebonato (1998) and Rebonato (1999) highlight 
the importance of jointly calibrating volatilities and coiTelations. These works, together with Longstaff, Santa-Clara & 
Schwartz (2000) also detail some of the most common non-convex calibration techniques based on parametrizations 
of the forward rates covariance factors. The mixed static-dynamic hedging formulation of the pricing problem has its 
source in the works by El Karoui, Jeanblanc-Picque & Shreve (1998), Avellaneda, Levy & Paras (1995) and Avellaneda 
& Paras (1996), while Romagnoli & Vargiolu (2000) provide some closed-form results in the multivariate case. 

Here, we show how the dual solution to the calibration program provides a complete description of the sensitivity 
to changes in market condition. In fact, because the algorithms used to solve the calibration problem jointly solve the 
problem and its dual, the sensitivity of the calibrated covariance matrix is readily available from the dual solution to 
the calibration program. When the objective in the calibration program is another swaption's price, the dual solution 
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also describes an approximate solution to the optimal hedging problem in AveUaneda & Paras (1996), which computes 
the price of a derivative product as the sum of a static hedging portfolio and a dynamic strategy hedging the worst-case 
residual risk. We also show how semidefinite programming can be used to efficiently solve the problem of optimally 
managing the Gamma exposure of a portfoho using vanilla options, as posed by Douady (1995). 

The results we obtain here underline the key advantages of applying semidefinite programming methods to the 
calibration problem: besides their numerical performance, they naturally provide some central results on sensitivity 
and risk-management. They can also eliminate numerical errors in sensitivity computations that were caused by the 
inherent instability of non-convex calibration techniques. 

The paper is organized as follows: In the next section, we quickly recall the calibration program construction 
for the Libor Market Model. Section three shows how to compute key sensitivities from the dual solution to this 
calibration problem. A fourth section details how these results can be used to form hedging portfolios. Finally, in the 
last section, we present some numerical results. 



2 Model Calibration 

In this section, we begin by briefly recalling the Libor Market Model construction along the lines of Brace et al. (1997) 
(see also Jamshidian (1997), Sandmann & Sondermaim (1997) and Miltersen, Sandmann & Sondermann (1997)). We 
then describe how to form the calibration program. 

2.1 Zero coupon dynamics 

We use the Musiela parametrization of the Heath, Jarrow & Morton (1992) setup. r{t, 9) is the continuously com- 
pounded instantaneous forward rate at time f , with duration 9. To avoid any confusion, Roman letters will be used for 
maturity dates and Greek ones for durations. The zero-coupon is here computed as 

B{t,T)=exp[- ' r{t,9)d9^ (1) 

All dynamics are described in a probability space (D,,{Ft;t > 0},P) where the filtration {Ft;t > 0} is the P- 
augmentation of the natural filtration generated by ad dimensional Brownian motion W = {Wt,t > 0}. The savings 
account is defined by: 

P, =exp^^ r{s,0)ds^ 

and represents the amount generated at time f > by continuously reinvesting 1 euro in the spot rate r{s, 0) during the 
period < i < As in Heath et al. (1992), the absence of arbitrage between all zero-coupons and the savings account 
states that: 

= B{0, T) exp j'^ a{s,T- s)dWs J^\a{s,T -s)\^ds^ (2) 

is a martingale under P for all T > 0, where for all > the zero-coupon bond volatility process {(7(f , 0); > 0} is 
fj-adapted with values in W^. We assume that the function 9 i — > (y{t, 9) is absolutely continuous and the derivative 
T{t, 9) = d/d9{(y{t, 9)) is bounded on x CI. 

2.2 Libor diffusion process 

All Heath et al. (1992) based arbitrage models are fully specified by their volatility structure and the forward rates 
curve today. The specification of the volatihty o{t, 9) in Brace et al. (1997) is based on the assumption that for a given 
underlying maturity 5 (for ex. 3 months) the associated forward Libor process {L{t,9y,t > 0} with maturity 9 defined 
by: 

l + 5L{t,9) = exp(^J ^ r{t,v)dvj 
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has a lognormal volatility structure. Using the Ito formula combined with the dynamics detailed above, we get as in 
Brace etal. (1997): 

dL{s,6) = (^^ + (2±il^o'^{s,e + 5){a\s,e + 5)-a\s,e)))ds 
+^ (1 + 5L{s,e)) (C7^(5, + 5) - C7«(5, e))dWs 

where Ws is the Brownian motion defined above and the deterministic volatiUty function y : R\ i — > M.^ is bounded and 
piecewise continuous. To get the desired lognormal volatility for Libors we must specify the zero-coupon volatility 

o(t,Ti) as: 

where {7i , . . . , 7m} is a calendar with period 5 and the Forward Rate Agreement K{t, T) = L{t, T — t) dynamics are 
given by: 

dK{s, T) = y{s, T - s)K{s, T) [a{s, T-s + 5)ds + dWs] 
Finally, as in Brace et al. (1997) we set a (f, 0) = for < < 5. 

2.3 Swaps 

A swap rate is the rate that zeroes the present value of a set of periodical exchanges of a fixed coupon against a floating 
coupon equal to a Libor rate. In a representation that is central in swaption pricing approximations, we can write swaps 
as baskets of forwards (see for ex. Rebonato (1998)). For example, in the case of a swap with start date Ts and end 
date Tiv: 

N 

swap{t, Ts, Tn) = £ (Oi{t)K{t, Tt) (5) 

(=5 

where Ts, Tf^ are calendar dates in {Ji , . . . , r^}: and: 

(di{t) = '^^siTi,Ti+i)B{t,Ti+i) 
Level (t, Ts,Tn) 

with cvg{Ti, Tj+i), the coverage (time interval) between Ti and 7^+i and Level(t,Ts,Ti^) the level payment, i.e. the sum 
of the discount factors for the fixed calendar of the swap weighted by their associated coverage: 

N 

Level{t,Ts,TN) = Y,cvg{Ti,Ti+i)B{t,Ti+i) 

i=S 

2.4 Swaption price approximation 

As in Brace & Womersley (2000), d'Aspremont (2003) or Kurbanmuradov et al. (2002), we approximate the swap 
dynamics by a one-dimensional lognormal process, assuming the weights (Oi{t) are constant equal to co;: 

dswap{s,Ts,T,) ^ f ^_.^(,^7._,y^.yL 
swap{s,T,TM) ^s 

where 

Kit,Ti) 

(Oi = (0, 



' swap{t,Ts,TN) 

is computed from the market data today and W,^^^ is a d dimensional Brownian motion under the swap martingale 
measure defined in Jamshidian (1997), which takes the level payment as a numeraire. We can use the order zero basket 
pricing approximation in Huynh (1994) and compute the price of a payer swaption starting with maturity T, written 
on swap{s,Ts,Ti\f), with strike K using the Black (1976) pricing formula: 

Levelit, Ts, T^) {swap{t, Ts, TN)N{h) -KN{h- V^^)) (8) 
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where 

h = — — — 

and swap{t,Ts, Jiv) is the value of the forward swap today with the variance Vt^ given by: 

(=S ' \iJ=S / 



^'Tr{Cl,rs)ds (9) 

We note S" the set of symmetric matrices of size nxn. This cumulative variance is a linear form on the forward rates 
covariance, with Q,t and e S^"'^"'"^ constructed such that: 

Cli = aoF = {6>i(i>j)ij^^s,N\ t and r, = {{y{sJi-s),y{sJj-s)))^ .^^^^^^ h 

where is the covariance matrix of the forward rates (or the Gram matrix of the y{s, Ti — s) volatiUty function defined 
above). Here swaptions are priced as basket options with constant coefficients. As detailed in Brace & Womersley 
(2000) or d' Aspremont (2003), this simple approximation creates a relative error on swaption prices of 1 — 2%, which 
is well within Bid- Ask spreads. Finally, we remark that caplets can be priced in the same way as one period swaptions. 

2.5 Semidefinite programming 

In this section we give a brief introduction to semidefinite programming. 
2.5.1 Complexity 

A standard form Unear program can be written: 

minimize c^x 
subject to Ax = b 

x^O 

in the variable x € M", where x ^ means here x componentwise nonnegative. Because their feasible set is the intersec- 
tion of an affine subspace with the convex cone of normegative vectors, the objective being Unear, these programs are 
convex. If the program is feasible, convexity guarantees the existence of a unique (up to degeneracy or unboundedness) 
optimal solution. 

The first method used to solve these programs in practice was the simplex method. This algorithm works well 
in most cases but is known to have an exponential worst case complexity. In practice, this means that convergence 
of the simplex method cannot be guaranteed. Since the work of Nemirovskii & Yudin (1979) and Karmarkar (1984) 
however, we know that these programs can be solved in polynomial time by interior point methods and most modem 
solver implement both techniques. 

More importantly for our purposes here, the interior point methods used to prove polynomial time solvability 
of Unear programs have been generalized to a larger class of convex problems. One of these extensions is caUed 
semidefinite programming. A standard form semidefinite program is written: 

minimize Tr(CX) 

subject to Tr(A,X) =bi, i = 1 , . . . ,m (10) 

in the variable X e S", where X^Q means here that X is positive semidefinite. Nesterov & Nemirovskii (1994) showed 
that these programs can be solved in polynomial time. A number of efficient solvers are available to solve them, the 
one used in this work is called SEDUMI by Sturm (1999). In practice, a program with « = 50 will be solved in less 
than a second. In what follows, we wiU also formulate semidefinite feasibility problems: 

find X 

subject to Tr(A,X) = bi, i = l,...,m 

X ^ 



4 



in the variable X gS". Their solution set is convex as the intersection of an affine subspace with the (convex) cone of 
positive semidefinite matrices and a particular solution can be found by choosing an objective matrix C and solving 
the corresponding semidefinite program (10). We will see below that most duahty results on hnear programming can 
be extended to semidefinite programs. 



2.5.2 Semidefinite duality 

We now very briefly summarize the duality theory for semidefinite programming. We refer again the reader to Nesterov 
& Nemirovskii (1994) or Boyd & Vandenberghe (2004) for a complete analysis. A standard form primal semidefinite 
program is written: 

minimize Tr{CX) 

s.t. Tr{AiX)=bi, i=l,...,m (11) 

XhO 

in the variable X e S''^. For X >r 0, y e W", we form the following Lagrangian: 

m 

L{X,y) = Tr(CX) + £ {h - Tr{AkX)) 
k=l 

(m \ m 

'£{-y,Ak+c)x] + '£ykh 
k=l J k=l 

and because the semidefinite cone is self-dual, we find that L{X,y) is bounded below in X ^ iff: 

m 

c-Y,ykAkhO 

k=\ 

hence the dual semidefinite problem becomes: 

maximize L^LiW^yt ^.j. 

When the program is feasible, most solvers produce both primal and dual solutions to this problem as well as a 
certificate of optimahty for the solution in the form of the associated duality gap: 

^i=Tr(x(c-pykA,}j^ 

which is an upper bound on the absolute error. If on the other hand the program is infeasible, the dual solution provides 
a Farkas type infeasibility certificate (see Boyd & Vandenberghe (2004) for details). This means that, for reasonably 
large problems, semidefinite programming solvers can be used as black boxes. 



2.5.3 Cone progranuning 

The algorithms used to solve linear and semidefinite programs can be extended a little further to include second order 
cone constraints. A cone program mixes linear, second order and semidefinite constraints and is written: 

minimize Tr(CX) 

subjectto \\d]vec{X)+ej\\<fJvec{X)+gj, j=l,...,p 
Tr{AiX)=bi, i=l,...,m 

xyo 

2 

in the variable X e S", where vec(X) turns X into a vector x e M" by stacking its columns. Again, a direct extension 
of the duality results above is valid for cone programs and solvers such as SEDUMI by Sturm (1999) give either both 
primal and dual solutions or a certificate of infeasibility in polynomial time. 
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2.6 The calibration program 

Here, we describe the practical implementation of the calibration program using the swaption pricing approximation 
detailed above. This is done by discretizing in s the covariance matrix F^. We suppose that the calibration data set 
is made of m swaptions with option maturity Tsy. written on swaps of maturity 7)v^ — Tg^. foxk= 1, ...,m, with market 
volatiUty given by a^. 

2.6.1 A simple example 

In the simple case where the volatihty of the forwards is of the form 7(s, T — s) = y{T — s) with y piecewise constant 
over intervals of length 5, and we are given the market price of k swaptions with a^T]^ G M.+ the Black (1976) cu- 
mulative variance of swaption k written on swap{t, Ts^. , Tf/^), the calibration problem becomes, using the approximate 
swaption variance formula in (9): 

find X 

s.t. Tr(%X) = a^Ti,, k=l,...,m (13) 
XhO 

which is a semidefinite feasibility problem in the covariance matrix X e S*' and Q.^ = Y^jLi ^fkj with (p/^j e S''^ the 
rank one matrix with submatrix d);td)J starting at element {j, j) and all other blocks equal to zero, with M = maxkNk. 
In other words 

with all other elements equal to zero. See Brace et al. (1999) or d'Aspremont (2003) for further details. 

2.6.2 The general case 

Here we show that for general volatiUties 7(5, 7" — i), the format of the calibration problem remains similar to that of 
the simple example above, except that X becomes block-diagonal. In the general non-stationary case where 7 is of 
the form 7(5, T — s) and piecewise constant on intervals of size 5, the expression of the market cumulative variance 
becomes, according to formula (9): 

(^kTs, = £ 5TV = / £ d>,d),- {Y{s,Ti-s),Y{s,Tj-s)) ds 
i=\ J' \ij=s ) 

where G S*'"' is a block-matrix with submatrix (h^ail starting at element [S^ — i,Sk — i) and all other blocks equal 
to zero. Here X,- is the Gram matrix of the vectors 7(7^, Tj — 7]). CaUbrating the model to the swaptions k= {1, ...,m) 
can then be written as the following semidefinite feasibihty problem: 

find Xi 

s.t. lfii5Tr(i2i,X) =0-275,, k=l,...,m 
XihO, i=l,...,M 

with variables X, G S'*^ ' . We can write this general problem in the same format used in the simple stationary case. Let 
X be the block matrix 



X = 



^1 


. 









. . 





. Xm 



the caUbration program can be written as in (13): 

find X 

s.t. TriClkX) = a^Ts,, k=l,...,m (14) 
X ^ 0, X block-diagonal 
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except that and X are here "block-diagonal" matrices. We can also replace the equality constraints with Bid- Ask 
spreads. In the simple case detailed in (13), the new calibration problem is then written as the following semidefinite 
feasibility problem: 

find X 

s-t. (^BidjcTs, < M^kX) < (yik,A, k = 1, ...,m 
X 1 

in the variable X e S*', with parameters £1^, <yBidv ^Askv ^h- ^S^^i' we can rewrite this program as a semidefinite 
feasibility problem: 

find X 

( \ Q-k 0^\ X Ol\ 
s.t. -n-l / f/i I =ol^,^Js, 



Tr 













X 











/ 
































f/2 _ 












" X 






















Ui 











-/ 













-'Bid,k 



Ts., k=l,...,m 



X,Ui,U2hO 

in the variables X, f/i , f/i € S*', which can be sunomarized as 

find X 

s.t. Tr{ClAsk,kX) = C72,i,i^i, 
'^{^Bid,kX) = Og^ji^ 

it yO, it block-diagonal 



Ts, k=l,.. 



(15) 



with it, Cli^ e S''*'^. Because of these transformations and to simplify the analysis, we will always discuss the stationary 
case with equality constraints (13) in the following section, knowing that all results can be directly extended to the 
general case using the block-diagonal formulation detailed above. 



3 Sensitivity analysis 

In this section, we show how the dual optimal solution can be exploited for computing solution sensitivities with 
minimiil numerical cost. 

3.1 Computing sensitivities 

Let us suppose that we have solved both the primal and the dual calibration problems above with market constraints 
a^Ts^ and let X°p' and y"?' be the optimal solutions. Suppose also that the market swaption price constraints are 
modified by a small amount u G M"*. The new caUbration problem becomes: 

minimize Tr{CX) 

s.t. Tr(akX) = a^Ts,+Uk k=l,...,m (16) 

X^O 

in the variable X gS'^ with parameters Q.ii,C gS'^ and o^Tsi^ + G Here C is, for example, an historical estimate 
of the covariance matrix. If we note />°p'(m) the optimal solution to the revised problem, we get (at least formally for 
now) the sensitivity of the solution to a change in market condition as: 

^=-,r (17) 

where y°P' is the optimal solution to the dual problem (see Boyd & Vandenberghe (2004) for details). As we will see 
in this section and the next one, this has various interpretations depending on the objective function. Here, we want to 
study the variation in the solution matrix X°P', given a small change u in the market conditions. 
We start with the following definition. 
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Notation 1 Let us suppose that we have solved the general calibration problem in (16), we call X°P' and y"?' the 
primal and dual solutions to the above problem with m = 0. We note 



2opt 



k=l J 



the dual solution matrix. As in Alizadeh, Haeberly & Overton ( 1998), we also define the symmetric Kronecker product 
as: 



1 



{P®Q)K:=- {PKQ^ + QKP^ 



We note A and A*, the linear operators defined by: 
A : S'^ — >W 

X^AX:=(Tr(n,X)),^i_, 



and its dual 



We now follow Todd & Yildirim (2001) to compute the impact AX on the solution of a small change in the market 
price data {uk)i^^i i-e. given u small enough we compute the next Newton step AX. Each solver implements one 
particular search direction to compute this step and we define a matrix M, with M = / in the case of the A.H.O. search 
direction based on the work by Alizadeh et al. (1998), see Todd & Yildirim (2001) for other examples. We also define 
the linear operators: 

£ = Z°P' ® M and F = MX"?' ® 7 



and their adjoints 



:Z°P'®MandX°P'M®/ 



We remark that if A,B e S commute, with eigenvalues a, j3 G M and common eigenvectors v, for i — 1, ...,M, then 
A®B has eigenvalues {atjij + aj^i) for ij = 1,...,M and eigenvectors v,vf if i = j and (vjvj + v^vf ) if ; ^ j for 
i,j = l,...,M. Provided the strict feasibility and nonsingularity conditions in §3 of Todd & Yildirim (2001) hold, we 
can compute the Newton step AX as: 



AX 



E-^FA* 



{AE-^FA 



, -1 



and this wiU lead to a feasible point X"?' + AX y iff the market variation movement u is such that: 



E-^FA* 



(ae-^fa* 



(X°P')" 



< 1 



(18) 



(19) 



where 1| • II2 is the I2 norm. The intuition behind this formula is that semidefinite programming solvers are based on the 
Newton method and condition (19) ensures that the solution X°P' remains in the region of quadratic convergence of the 
Newton algorithm. This means that only one Newton step is required to produce the new optimal solution X^P' + AX 
and (18) simply computes this step. The matrix in (18) produces a direct method for updating X which we can now 
use to compute price sensitivities for any given portfolio. 

This illustrates how a semidefinite programming based calibration allows to test various realistic scenarios at a 
minimum numerical cost and improves on the classical non-convex methods that either had to "bump the market data 
and recalibrate" the model for every scenario with the risk of jumping from one local optimum to the next, or simulate 
unrealistic market movements by directly adjusting the covariance matrix. One key question is stability: the calibration 
program in (13) has a unique solution, but as usual this optimum can be very unstable and the matrix in (18) badly 
conditioned. In the spirit of the work by Cont (2001) on volatility surfaces, we look in the next section for a way to 
solve this conditioning issue and stabilize the calibration result. 



3.2 Robustness 

The previous sections were focused on how to compute the impact of a change in market conditions. Here we will focus 
on how to anticipate those variations and make the calibrated matrix robust to a given set of scenarios. Depending on 
the way the perturbations are modelled, this problem can remain convex and be solved very efficiently. Let us suppose 
here that we want to solve the caUbration problem on a set of market Bid- Ask spreads data. A direct way to address 
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the conditioning issues detailed in the last section is to use a Tikhonov stabilization of the caUbration program as in 
Cont (2001) and solve the following cone program: 



minimize 
s.t. 



Il^l 

r2 



<d,kTs, < Tr(n,X) < al,^,Ts„ k = 1, 



(20) 



in the variable X e S'^ with parameters i2i,C € S'^ and (ysidk'^Si,, ^Askk^h ^ ^® absence of any information 

on the uncertainty in the market data, we can simply maximize the distance between the solution and the market 
bounds to ensure that it remains valid in the event of a small change in the market variance input. As the robustness 
objective is equivalent to a distance maximization between the solution and the constraints (or Chebyshev centering), 
the input of assumptions on the movement structure is equivalent to a choice of norm. Without any particular structural 
information on the volatility market dynamics, we can use the /<» norm and the calibration problem becomes: 



maximize t 

Ts,+t<Tr{akX)<al,^^Ts,- 
X>0 



-t, k = l,...,m 



in the variables ? e IR and X & S . Using the li norm instead, this becomes: 



maximize 
s.t. 



^m,kTs, + tk < Tr(n^X) < crj -f*, k = l,...,m 



xyo 

in the variables ti,...,tm € M and X G S''^. The problems above optimally center the solution within the Bid- Ask 
spreads, which makes it robust to a change in market conditions given no particular information on the nature of that 
change. In the same vein, Ben-Tal, El Ghaoui & Lebret (1998) also show how to design a program that is robust to a 
change in the matrices Q.^. However, because the matrices Q.^ are computed from ratios of zero-coupon bonds, their 
variance is negligible in practice compared to that of a^. 

Suppose now that V is a statistical estimate of the daily co variance of c^Tsj. (the mid-market volatilities in this 
case) and let us assume that these volatihties have a Gaussian distribution. We adapt the method used by Lobo, 
Vandenberghe, Boyd & Lebret (1998) for robust L.P. We suppose that the matrix V has full rank. We can then 
caUbrate the model to this information: 



maxrmize 
s.t. 



Tr{CX) 
Tri^kX) 

xyo 



<<i) 1 



forA: = 1, ...,m 



where || • ||^ is the Zoo norm and <I>(x) is given by 

1 f-^ 

<^{x) = 1 = / exp{-u^/2)du 

There is no guarantee that this program is feasible and we can solve instead for the best confidence level by forming 
the following program: 

minimize t 

s.t. Tr (iliX) - a^Ts, =Vk foTk= l,...,m 



xyo 



<t 



The optimal confidence level is then tj = ^{t) and "centers" the calibrated matrix with respect to the uncertainty in 
o^Tsj^. This is a symmetric cone program and can be solved efficiently. 
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4 Hedging 



In this section, we show how semidefinite programming calibration techniques can be used to build a superreplicating 
portfolio, approximating the upper and lower hedging prices defined in El Karoui & Quenez (1995). An efficient tech- 
nique for computing those price bounds with general non-convex payoffs on a single asset with univariate dynamics 
was introduced in AveUaneda et al. (1995) and recent work on this topic by RomagnoU & Vargiolu (2000) provided 
closed-form solutions for the prices of exchange options and options on the geometrical mean of two assets. 



4.1 Approximate solution 

Here, using the approximation in (8), we first compute arbitrage bounds on the price of a basket by adapting the 
method used by AveUaneda et al. (1995) in the one-dimensional case. We then provide approximate solutions for 
these arbitrage bounds on swaptions and show how one can use the dual solution to build an optimal hedging portfoUo 
in the sense of AveUaneda & Paras (1996), using derivative securities taken from the calibration set. 

As in AveUaneda & Paras (1996), the price here is derived from a mixed static-dynamic representation: 

Price = Min {Value of static hedge -|- Max (PV of residual liability) } (21) 

where the static hedge is a portfoUo composed of the calibration assets and the maximum residual liability is computed 
as in El Karoui et al. (1998) or AveUaneda et al. (1995). This hedging representation translates into the pricing problem 
the market habit of calibrating a model on the instrument set that will be used in hedging. The static portfoUo uses 
these instrument to reduce the payoff risk, while the dynamic hedging part hedges the remaining risk in conservative 
manner. 

Furthermore, because of the sub-additivity of the above program with respect to payoffs, we expect this diversifi- 
cation of the volatiUty risk to bring down the total cost of hedging. Let K{t) = {K{t, Ti))^^^ ^ and suppose we have 
a set of market prices Q for k— 1 , . . . ,m, of swaptions with corresponding market volatilities a^, coefficient matrices 
Slk S S*' and payoffs hi^{K{t)). As in AveUaneda & Paras (1996) we can write the price (21) of an an additional 
swaption with payoff ho{K{t)): 



inf \ y AfcC;t+ supE 



k=i 



(22) 



where P varies within the set of equivalent martingale measure and j3 (T) is the value of the savings account in T. We 
can rewrite the above problem as: 

inf jsup (^E'' [l3{T)-%iK{T))] - £ A, [l5iT)-'h,{K{T))] -Q) j | 

where we recognize the optimum hedging portfoUo problem as the dual of the maximum price problem above: 

maximize £^ 
s.t. E'' 



P{T)-^ho{K{T)) 



= Ck, k=l,...,m 

Using (8), we get an approximate solution by solving the following problem: 

maximize BSo{Tr{Q.oX)) 

s.t. BSk{Tr{Q.kX))=Ck, k=l,...,m 

and its dual: 

inf I sup (^BS{Tr{^oX)) -f^h (BS Cn-(il^Z)) - Q) j | 

The primal problem, after we write it in terms of variance, becomes the foUowing semidefinite program: 

maximize a^J' = 'Tr{CloX) 
s.t. Tr(D.kX) = a^Tk, k=l,...,m 

X^O 
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Again, we note y°^^ G W" the solution to the dual of this last problem: 

minimize Y!k=iyk<^k'^k 



The KKT optimality conditions on the primal-dual semidefinite program pair above (see Boyd & Vandenberghe (2004) 
for example) can be written: 

= l.'^^iyk^kX-aoX 
Tr{akX) = G^Tk, k=l,...,m 
0-<X 

and we can compare those to the KKT conditions for the price maximization problem: 

\ xz = o " 

] B5i(Tr(nfcX))=q, k=l,...,m 

[ o^x,z 

with dual variables A G K™ and Z G S„. An optimal dual solution k^^^ for the price maximization problem can then be 
constructed from y°P', the optimal dual solution of the semidefinite program on the variance, as: 

, opt ^ _ opt dBSo{Tr{^oX))/dv 
^ dBSk{Tr{akX))/dv 

which gives the amount of basket k in the optimal static hedging portfoho defined by (2 1) . 



4.2 The exact problem 

The bounds found in the section above are only approximate solutions to the superreplicating problem. Although the 
relative error in the swaption price approximation is known to be about 1-2%, it is interesting to notice that while being 
somewhat intractable, the exact problem shares the same optimization structure as the approximate one. Let us recall 
the results in Romagnoli & Vargiolu (2000). If, as above, we note C{K{t)^t) the superrepUcating price of a basket 
option, then C{K{t),t) is the solution to a multidimensional Black-Scholes-Barenblatt (BSB) equation. We can create 
a superrephcating strategy by dynamically trading in a portfoho composed of A} = ^{t,K{t,Ti)) in each asset. The 
BSB equation in Romagnoli & Vargiolu (2000) can be rewritten in a format that is similar to that of the approximate 
problem above, to become: 

dC(x,t) I „^ Tv/'^- ( \ d'^C(x.t) 



dt 



2 maxrgATr {diag{x) ^-^^diag{x)Tj = 



[ C{x,T) = {r:^,(OiXi-k)+ 

where diag{x) is the diagonal matrix formed with the components of x and F = 77^ is the model covariance matrix. 
If the set A is given by the intersection of the semidefinite cone (the covariance matrix has to be positive semidefinite) 
with a polyhedron (for example approximate price constraints, sign constraints or bounds on the matrix coefficients, 
...), then the embedded optimization problem in (Romagnoli & Vargiolu (2000)) becomes a semidefinite program: 

f^, , ,d^C(x,t) , , , 
miucTr \ rdiag(x) — ^-^ — diag{x) 

on the feasible set A. We recover the same optimization problem as in the approximate solution found in the section 

above, the only difference being here that the solution to the general problem might not be equal to a Black-Scholes 
price. This gives a simple interpretation of the embedded optimization problem in the BSB equation developed in 
Romagnoli & Vargiolu (2000). 



11 



4.3 Optimal Gamma Hedging 

For simplicity here, we work in a pure equity framework and, along the lines of Douady (1995), we study the problem 
of optimally adjusting the Gamma of a portfolio using only options on single assets. This problem is essentially 
motivated by a difference in liquidity between the vanilla and basket option markets, which makes it impractical to 
use some baskets to adjust the Gamma of a portfolio. Suppose we have an initial portfolio with a Gamma sensitivity 
matrix given by F in a market with underlying assets xi for j = 1, . . .n. We want to hedge (imperfectly) this position 
with yi vanilla options on each single asset Xi with Gamma given by Yi. We assume that the portfolio is maintained 
delta-neutral, hence a small perturbation of the stock price will induce a change in the portfoUo price given by: 

AP(x + Ax) = P{x) + i A.v^r(y)Ax 

where r(y) = r + diag{Y)y, with diag(7) the diagonal matrix with components Yi. As in Douady (1995), our objective 
is to minimize in y the maximum possible perturbation given by: 

max |Ax^r(y)Ax| 

where is the elhpsoid defined by 

S'={ueR"\u'^I.u = l} 

with Zy = cov{xi,Xj) for i,j = l,...,n, the covariance matrix of the underlying assets. This amounts to minimizing 
the maximum eigenvalue of the matrix Tr{y) and can be solved by the following semidefinite program: 

minimize t 

subject to -// ^ Zr + Idiag{Y)diagiy) ^ tl 

in the variables f G R and j e R". 

5 Numerical results 

5.1 Price bounds 

We use a data set from Nov. 6 2000 and we plot in figure (1) the upper and lower bounds obtained by maximizing (resp. 
minimizing ) the volatility of a given swaption provided that the Libor covariance matrix remains positive semidefinite 
and that it matches the calibration data. We calibrate by fitting all caplets up to 20 years plus the following set of 
swaptions: 2Y into 5Y, 5Y into 5Y, 5Y into 2Y, lOY into 5Y, 7Y into 5Y, lOY into 2Y, lOY into 7Y, 2Y into 2Y. This 
choice of swaptions was motivated by liquidity (where all swaptions on underlying and maturity in 2Y, 5Y, 7Y, lOY 
are meant to be Uquid). For simpUcity, all frequencies are annual. For each [Maturity, Underlying) pair in figure (1), 
we solve the semidefinite program detailed in (13): 

find X 

s.t. Tr(%X) = olTs,, k = 1, ...,m 
X^O 

in the variable X e S^^, where (y^ and are computed as in (13) using the data in Table (1) and (2). The programs are 
solved using the SEDUMl code by Sturm (1999). We then compare these upper and lower bounds (dotted lines) with 
the actual market volatihty (solid line). Quite surprisingly considering the simphcity of the model (stationarity of the 
sliding Libor dynamics L(f, 0)), figure (1) shows that all swaptions seem to fit reasonably well in the bounds imposed 
by the model, except for the 7Y and lOY underlying. This is in line with the findings of Longstaff et al. (2000). 

5.2 Super-replication & calibration stability 

Here, we compare the performance and stabihty of the various calibration methodologies detailed here and in Rebonato 
(1999). For simplicity, we neglect the different changes of measure between forward measures and work with a 
multivariate lognormal model where the underlying assets xif for j = 1 , . . . , 5 follow: 

1 r 

dxi^s/xi^s = OidWs - -C7; Oids (23) 
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Figure 1: Calibration result and price bounds on a "Sydney opera house" set of swaptions. 



CapletVols(%,lYtolOY) 14.3 15.6 15.4 15.1 14.8 14.5 14.2 14.0 13.9 13.3 
CapletVols(%,llYto20Y) 13.0 12.7 12.4 12.2 12.0 11.9 11.8 11.8 11.7 12.0 



Table 1: Caplet volatilities. 



Swaption 


Vol (%) 


Cbi 










2Y into 5Y 


12.4 


0.22 


0.20 


0.20 


0.19 


0.18 


5Y into 5Y 


11.7 


0.22 


0.21 


0.20 


0.19 


0.18 


5Y into 2Y 


14.0 


0.51 


0.49 








lOY into 5 Y 


10.0 


0.22 


0.21 


0.20 


0.19 


0.18 


7Y into 5Y 


11.0 


0.23 


0.21 


0.20 


0.19 


0.18 


lOY into 2Y 


12.2 


0.51 


0.49 








lOY into 7Y 


9.6 


0.17 


0.16 


0.15 


0.14 


0.13 


2Y into 2Y 


14.8 


0.52 


0.48 









Table 2: Swaption volatilities and weights (data courtesy of BNP Paribas, London). 
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where Wg is a B.M. of dimension 5. At time 0, we set Xifi = .1 and the covariance matrix a^a is given by: 

/ 0.0144 0.0133 0.0074 0.0029 0.0013 \ 

0.0133 0.0225 0.0151 0.0063 0.0032 

0.0074 0.0151 0.0144 0.0065 0.0032 

0.0029 0.0063 0.0065 0.0081 0.0027 

\ 0.0013 0.0032 0.0032 0.0027 0.0036 / 

We then compare the performance of a delta hedging strategy implemented using various calibration techniques. In 
each case, we hedge a short position in an ATM basket option with coefficients wq = (0,0, .4, .1, .5), cahbrating the 
model on all single asset ATM calls and an ATM basket with weight ws = (1 , 1 , 0, 0, 0), the weights wi , . . . , W5 being 
then equal to the Euclidean basis. All options have maturity one year and we rebalance the hedging portfolio 33 times 
over this period. At each time step, we calibrate to these option prices computed using the model in (23). To test the 
stability of the caUbration techniques, we add a uniformly distributed noise to the calibration prices with amphtude 
equal to ±10% of the original price. 

We use four different cahbration techniques to get the covariance matrices and compute the deltas. In the first 
one, we use the exact model covariance above to compute baseUne results. In the second one, we use the simple 
stabilization technique detailed in (20) and solve: 

minimize 11.^112 

s.t. Tr{£lkX) = a^Ts,, k=l,...,6 (24) 

in the variable X, where and the matrices for ^ = 1, ... ,6 are computed from wi,. . . ,W6 as in (13). In the third 
one, we use the super replication technique detailed in §4.1 and calibrate the covariance matrix by solving 

maximize Tr(i2o^) 

s.t. TriClkX) = a^Ts,, k=l,...,6 (25) 

x-^o 

in the variable X, where and the matrices £2^ for A: = 0, ... ,6 are computed from wq, . . . as in (13). Finally, 
our fourth caUbrated covariance matrix is caUbrated using the two factor parametrized best fit technique detailed in 
Rebonato (1999). 

For each calibration technique, we record the ratio of the delta hedging portfoho's P&L to the initial option pre- 
mium at every time step. In figure (2), we plot the P&L distributions for the super-hedging strategy (25) and the 
Rebonato (1999) calibration technique. We notice that while discrete hedging error makes the super-replication some- 
what imperfect, the super-hedging calibration technique has a much higher P&L on average than the best fit cahbration 
in Rebonato (1999). Furthermore, the super-hedging calibration produces a positive P&L in 68% of the sample sce- 
narios, while the parametrized calibration has a positive P&L in only 41% of them. 

Finally, in table (2) we detail some summary statistics on the relative P&L of the four cahbration techniques 
detailed above and on the change in covariance matrix between time steps. In this table, "mean change" is the average 
norm of the change in the calibration matrix at each time step, while "Mean P&L" and "StDev" are the mean and 
standard deviation of the ratio of the hedging portfolio's P&L to the option premium. Finally, the "real covariance" 
calibration technique uses the actual model covariance matrix to compute the delta, the "robust" technique uses the 
solution to (24), the "super-hedging" one uses the solution to (25), while the "parametrized" uses the parametric 
algorithm described in Rebonato (1999). 

Calibration method Real covariance Robust Super-Hedging Parameterized 
Mean P&L ^0002 0083 Ol37 ^0109 
StDev P&L 0.253 0.338 0.344 0.316 

Mean change (10-5) 3.11 3.36 4.77 

Table 3: Hedging P&L and covariance stability statistics for various calibration techniques. 

We notice that, as expected, the super hedging strategy improves the average P&L. However, while the robust 
calibration algorithm does show a smaller average change in calibrated covariance matrix, this does not translate 
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Figure 2: P&L distribution for the super-hedging (light) and the parametric (dark) calibration techniques. 

into significantly smaller hedging P&L standard deviation. This is perhaps due to the fact that the basket options 
considered here are not sensitive enough for these changes in the covariance to have a significant impact on the 
hedging performance. 

6 Conclusion 

The results above have showed how semidefinite programming based calibration methods provide integrated cali- 
bration and risk-management results with guaranteed numerical performance, the dual program having a very natural 
interpretation in terms of hedging instruments and sensitivity. Furthermore, these techniques make possible the numer- 
ical computation of super-hedging strategies (detailed in section 4.1), which seem to perform well in our simulation 
examples. 
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